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Abstract 



We present results concerning the non-perturbative evaluation of the renormalisation 
constant for the quark field, Zg, from lattice simulations with twisted mass quarks and three 
values of the lattice spacing. We use the RI'-MOM scheme. Zq has very large lattice spacing 
artefacts ; it is considered here as a test bed to elaborate accurate methods which will be used 
for other renormalisation constants. We recall and develop the non-perturbative correction 
methods and propose tools to test the quality of the correction. These tests are also applied 
to the perturbative correction method. We check that the lattice spacing artefacts scale 
indeed as a^p^. 

We then study the running of Zq with particular attention to the non-perturbative effects, 
presumably dominated by the dimension- two gluon condensate (^^) in Landau gauge. We 
show indeed that this effect is present, and not small. We check its scaling in physical units 
confirming that it is a continuum effect. It gives a ~ 4% contribution at 2 GeV. 

Different variants are used in order to test the reliability of our result and estimate the 
systematic uncertainties. 

Finally combining all our results and using the known Wilson coefficient of (A^) we find 
5^(/i^)(j4^)^2 CM = 2.01(11) (lo;73) GeV^ at /x = 10 GeV, in fair agreement within uncertain- 
ties with the value indepently extracted from the strong coupling constant. We convert the 
non-perturbative part of Zq from RI'-MOM to MS. Our result for the quark field renormalisa- 



tion constant in the MS scheme is p'''*((2 GeV)^, c/^^^J = 0.750(3)(7) - 0.313(20) (^barc" 
1.5) for the perturbative contribution and ^J5Snon-pcrturbativc^^2 GeV)^, ^^^^J = 0.781(6)(21)- 
0.313(20) {gl^^^ - 1.5) when the non-perturbative contribution is included. 
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1 Introduction 

Computing matrix elements in lattice Quantum ChromoDynamics (LQCD) needs often the com- 
putation of renormalisation constants. Indeed, even if the lattice computation contains only O(a^) 
lattice artefacts, the bare quantities differ from the continuum ones by 0{g'^) ~ 0(1/ log(a^)) which 
is unacceptable. Renormalisation restores the O(a^) accuracy. It is also known since long that 
these renormalisation constants need to be computed non-perturbatively, using LQCD techniques. 
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Several non-perturbative methods have been proposed. Let us here concentrate on those based 
on the MOM scheme. They start from the computation of Green functions of quarks, gluons and 
ghosts at large enough momenta in a fixed gauge, usually the Landau gauge. This gives the 
renormalisation constant Z{fi) at many values of the scale /i. Assuming that our goal is to deliver 
the renormalisation constant in the MS scheme at say 2 GeV (a typical phenomenological scale), 
one must then use results from perturbative QCD to convert MOM into MS and run to 2 GeV. 
The running of ZMOAf(/^) is a very powerful testing tool: indeed perturbative QCD is only useful if 
we are in the perturbative regime, i.e. at large enough momenta. The only way to check whether 
this is the case is to compare lattice data with the perturbative running. In this framework, it 
turns out that this is not always the case. 

Deviations from perturbative running can be analysed via Wilson operator expansion and the 
Shifman-Vainshtein-Zakharov (SVZ) sum rules. It turns out that the dominant non-perturbative 
correction in Landau gauge is due to the non vanishing vacuum expectation value of the only 
dimension-two operator : = A'^^A"'^ and that it is not small [2]-[8j. It is thus necessary 
to carefully look for the possibility of such a contribution, which appears in the OPE, as a 
contribution up to logarithmic corrections. The coefficient of this contribution is equal 
to the vacuum expectation value g'^{A^) times a Wilson coefficient that has to be computed in 
perturbation theory, and has been up to three loops for propagators [o] . To argue that a measured 
l/p^ contribution is a continuum power correction and not a lattice artefact, we must check that the 
term in the fit scales with lattice spacing when expressed in physical units. To further argue 
that this is indeed due to {A^) we must compare the resulting {A^) from different quantities and 
thus check the universality of the condensates which is on the ground of the SVZ technology. The 
theory of Wilson operator expansion is then constraining: since there exists only one dimension- 
two operator, provided that it is renormalized with the same prescription for all these quantities, 
all the different estimates of {A^) should coincide within errors, up to corrections, for a given 
value of Nf and of the dynamical masses. Of course, its extraction needs the coefficients of the 
Wilson expansion which are computable in perturbation. To test this universality of the extracted 
{A^) is one of the goals of our program of analysing many different quark and gluon quantities 
obtained from lattice gauge configurations produced by the European twisted mass collaboration 



(ETMC). We have also applied a criterium proposed in 10 to validate the way we use operator 
expansion. This paper makes one of the first steps in such a program. 

It is worthwile also to mention that several authors have elaborated further on the relation 
between this gauge-dependent gluon condensate, obtained in the Landau gauge, and possible 
1/p^-terms in gauge invariant quantities, and thus on the phenomelogical implications, mainly in 



connection with confinement scenarios, of such a dimension-two condensate 11 . 

All this can only be done once the lattice artefacts are eliminated or at least under control. The 
O(a^) artefacts can be quite large since we consider large momenta, while finite volume artefacts 
are minor. There are two main types of 0{a?) artefacts: 0{a?p^) artefacts which respect the 
continuum 0(4) rotation symmetry, and hypercubic artefacts which respect the hypercubic 
symmetry group but not 0(4). The latter are effects of the hypercubic symmetry of the lattice 
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action. We will identify the 0{a?p^) artefacts non-perturbatively by doing a fit of the running Z{^) 
which will include the perturbative running, the {A^) power correction and a term proportional 
to a^p^. Notice that, while the perturbative and {A^) running contributions must approximately 
scale in physical units, the artefacts must scale in lattice units. This is an additional check we 
shall perform. 

Concerning the elimination of hypercubic artefacts, which is better done before the above 
mentioned running fit, several methods have been proposed in literature: the democratic selection, 
the perturbative correction and the non-perturbative "egalitarian" one ( "egalitarian" because all 
the points are used on the same footing in this approach). We will discuss this in some detail later 
and perform extensive comparisons. In particular we will use a new quality test which consists 
in watching to what extent the "half-fishbone" structure, which raw lattice results for Zq always 
exhibit and which is a dramatic illustration of hypercubic artefacts, are corrected by every method. 

Although all the issues raised here concern all the renormalisation constants as well as the 
QCD coupling constant, we will concentrate in the following on Zq, that renormalises the quark 
field 



ZI'\b (1.1) 



where qs {qii) is the bare (renormalised) quark field. In the RI'-MOM scheme Zq is defined by 

Ul^' = V") = ^.'^riS.iM f] (1-2) 

where S'bare(p) is the bare quark propagator. Our goal is to compute that constant from LQCD 
with twisted Wilson quarks. 

In |4)|l2| astudy^of was performed from LQCD in the case Nj = using both the overlap 



and Wilson clover fermions. In 12 the exceptional size of hypercubic artefatcs was stressed and a 
non-perturbative elimination of hypercubic artefacts performed along the same principle as what 
is used here. In ji] the Wilson coefficient of the {A'^) was computed up to the leading logarithm 
approximation and applied to estimate the condensate from the LQCD data. The outcome was 
that a significant non-perturbative contribution from {A"^) was needed to account for the results. 
Notice that we do not expect the {A"^) to be similar or even close in the cases of Nf = and 
Nf = 2. 

Summarising the above discussion, we do here concentrate on Zq because we consider it as a 
kind of benchmark for the following reasons: 

• It has specially large hypercubic artefacts and is thus a good test bed for a correct treatment 
of these. 

• It has a vanishing anomalous dimension at leading order in the Landau gauge: it's pertur- 
bative running is thus soft. 



^In [ij Zq was denoted 
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• The Wilson coefficient of the (A^) condensate is rather large [4||9], which is an incentive to 
look carefully for non-perturbative contributions. 

In this paper, in order to test deeply the reliability of our results, we will compare many fits: 
perturbative/non-perturbative hypercubic correction, one- window/sliding- windows non perturba- 
tive hypercubic correction, effect of the total fitting range, fitting separately every /3 and global 
fit. As a consequence we will proceed as follows: 

• We will recall some general formulae concerning the perturbative and non-perturbative run- 
ning in the continuum 

• describe our lattice setting and our non-perturbative "egalitarian" method to eliminate 
hypercubic artefacts; 

• present the results concerning the perturbative correction method for hypercubic artefacts 
and show the quality checks; 

• present the results concerning the non-perturbative method to correct for hypercubic artefacts 
and show the quality checks, propose two types of fits, the sliding windows fit (SWF) and 
the one- window fit (OWF). 

• We will perform the running fit on the output of all the previously mentioned hypercubic 
corrected data, compare the results for the g'^{A^) for all these fits, check the scaling of 

• check the scaling of the oc a^p^ artefacts; 

• check the lattice spacing dependence of Z^"^^, the perturbative contribution to Zq, (x g"^; 

• study the range of variation of our results for g'^{A^) from the egalitarian method with 
one/sliding window(s), with varying fitting ranges, and the perturbative method with two 
realisations. We extract from there the systematic uncertainty. 

• We will join in one plot the three /3's and perform the fit of the running, 

• compare the resulting g'^{A^) with the one extracted from the strong coupling constant 
and with quenched estimates, and tested our procedure according to Martinelli-Sachrajda's 
criterium 
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conclude. 
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2 Running of 



2.1 Perturbative running 

In Landau gauge Zg has a vanishing anomalous dimension to leading order, i.e. its running starts 



at O(a^). The perturbative running has been computed up to four loops 13 and references 



therein. The needed formulae are accessible on the web site also indicated in ref. 13 



2.2 Wilson expansion and non-perturbative running 



To handle non-perturbative corrections we use the theory of Operator Product Expansion 14 and 



its application in estimating power suppressed non-perturbative corrections via vacuum expecta- 



tion values 15 . In Landau gauge there exists only one dimension-2 operator allowed to have a 



vacuum expectation value: = YllJil ^u^"^- The Wilson coefficient of this operator has been 



computed to leading logarithm in |4| and extensively for all propagators up to 0{a ) in [9|. 
2.2.1 (A^) tree level Wilson coefRcients for Zq 

In order to give a hint let us just sketch the tree level calculation of that Wilson coefficient. 

P P P 




Consider the above diagram, describing a quark propagating in a constant background gauge 
field. As a consequence the red bubble represents the interaction of the quark with this background 
field: igXa/2 jjL"". The Feynman rules are then applied as usual. Neglecting the quark mass it gives 

^ \l_i=l,a=l ^t'=l,a'=l ^ / 

where ( ) represents the vacuum expectation value, = Cp = 4/3 (proportional to the 

identity matrix in color space), the sum over fi gives a factor 4, and 

((A^)^) = (AV32), ((A-p)2) = (AV4) (2.2) 

from the homogeneity of the vacuum for rotations in space-time and color space. 

For Zg defined by Eq. (1.2), we get at tree level the following non-perturbative contribution 
due to {A^): 

12 p 
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2.2.2 The Wilson coefficients at 0{a^) 

The Wilson coefficient of (A^) for the quark propagator has been computed up to 0{a'^) in [o] in 
the MS scheme. Our lattice data refer naturally to the RI'-MOM scheme. Some work is needed 
to derive the correct analytic formula which allows a fitting of our lattice data. We have derived 
this in the appendix |Aj 



3 The lattice computations 

The results presented here are based on the gauge field configurations generated by the European 
Twisted Mass Collaboration (ETMC) with the tree-level improved Symanzik gauge action fT6] 
and the twisted mass fermionic action [T7| at maximal twist. 



3.1 The lattice action 

A very detailed discussion about the twisted mass and tree-level improved Symanzik gauge actions, 
and about the way they are implemented by ETMC, can be found in refs. 18 -21 . Here, for the 
sake of completeness, we will present a brief reminder of the twisted action and the run parameters 
for the gauge configurations that will be exploited in the present work (See tab. [T|. 

The Wilson twisted mass fermionic lattice action for two flavours of mass degenerate quarks, 
reads (in the so called twisted basis fT7l|22] ) 



(3.1) 
nv 

= -7. (V, + v;) - y v,v; , 

where mo is the bare untwisted quark mass and /ig the bare twisted quark mass, is the third 
Pauli matrix acting in flavour space and r is the Wilson parameter, which is set to r = 1 in the 
simulations. The twisted Dirac operator is defined as 

Aw = D^^ + mo + ij^TsHg . (3.2) 

The operators and V* stand for the gauge covariant nearest neighbour forward and backward 
lattice derivatives: 

D, = ^ [V,{x, y) + V;:(x, y)] = \ [5,,,.+^ W[x, y) - 6y,,.f, W{x, y)] ; (3.3) 

defining the operator as the discretized covariant derivative. The bare quark mass mo is related 
as usual to the so-called hopping parameter k, hj k = l/{8 + 2amo). Twisted mass fermions are 
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said to be at maximal twist if the bare untwisted mass is tuned to its critical value, merit- This is 
in practice done by setting the so-called untwisted PCAC mass to zero. 

In the gauge sector the tree-level Symanzik improved gauge action (tlSym) [16] is applied. 
This action includes besides the plaquette term U^^^^^^ also rectangular (1 x 2) Wilson loops U^^^^^- 
It reads 

= f E ( ^0 E {1 - KeTr(f/i^^;j}+6i ^ {1 - MeTr(f/2f J} j , (3.4) 

X \ ii,v = l fJ,,U = l / 

where (3 = G/g^, go being the bare lattice coupling and it is set 6i = —1/12 (with 6o = 1 — 8bi as 
dictated by the requirement of continuum limit normalization). Note that at 6i = this action 
becomes the usual Wilson plaquette gauge action. The run parameters for (3 and fig of the gauge 
configurations that will be exploited in the following can be found in tab. [T] 



/3 


a fm 


a-i GeV 




Volume 


# confs 


3.9 


0.083 


2.373 


0.004 


243 X 48 


100 


4.05 


0.0675 


2.897 


0.006 


243 X 48 


100 


4.2 


0.055 


3.58 


0.002 


243 X 48 


100 



Table 1: Run parameters of the exploited data from ETMC collaboration for the present study of 
Zg. The second column lists the lattice spacings which we have used in this study. It is easy to 
convert it to other lattice spacings. 



3.2 The computation of the quark propagator 



Computing the renormalisation constants for the quark propagator and the operators containing 
quark fields demands to compute first the gauge-fixed 2-point quark Green functions from the 
lattice. We exploited ETMC gauge configurations [23] obtained for /3 = 3.9, (3 = 4.05 and (3 = 4.2. 
After checking the small dependence of Zg on the dynamical and valence quark masses we decided 
to use only one mass for every (3, table [1} The lattice gauge configurations are transformed to 



Landau gauge by minimising the following functional of the SU(3) matrices, f/^(x) 



Fu[g] = Re 



X fj, ^ 



(3.5) 



with respect to the gauge transform g, by applying a combination of overrelaxation algorithm and 
Fourier acceleration |3 



^We end when < 10 and when the spatial integral of Aq is constant in time to better than 10 
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We compute quark propagators with a local source taken at a random point Xq on the lattice, 
in order to reduce the correlation between successive configurations: 

5(l/,Xo)r^'"'^" = D;^\y,xr^-^''^-'^'^so'/ix,Xo) , sof (x,Xo) = 5.,.o'^Mo^/3,/3o ; (3.6) 

where the equation is solved for every Bq = 1,3 and /3o = 1,4, and j = u,d labels the isospin. We 
perform the Fourier transform which is a 12 x 12 complex matrix 

5,(p)^5^e-^^(^-^'°)^,(l/,xo) . (3.7) 
y 

This is the Fourier transform of the quark incoming to the source (the arrow pointing towards the 
source). The Fourier transforms of the quark outgoing from the source is 

Sr{p)=l,SUph,, (3.8) 



where u ^ d]d = u. From Eq. (1.2) the lattice quark renormalisation constant Zq is given by 



where < ... > means here the average over the chosen ensemble of thermalised configurations and 
= - sinap^. The reason to use = - sinap^ is to get = 1 for a free fermion, or in other 
words, to ehminate hypercubic artefacts at tree level. 

3.3 The method of non-perturbative Hypercubic H{4:) correction 

The lattice estimates of the quark field renormalisation constant and the vertex functions lead to 
dimensionless quantities that, because of general dimensional arguments, depend on the strong in- 
teraction scale Aqcd and on the lattice momentum ap^. We have computed the Fourier transforms 
for the following momenta: 

P. = ^ n, = -iV,/4,...,iV,/4, P4 = ^^^^^;^ n, = -NT/A,...,Nr/A; 

(3.10) 

where pi = 1,3 are the spatial momenta and the time like. The antiperiodic boundary condition 
in the time direction explains the 7r(2n4 + 1) factor. 



The lattice action Eq. (3.1) and Eq. (3.4) is invariant under the hypercubic group H{A). 
However the boundary conditions and the difference between the spatial size Nl and the time-like 
one Nt = 2Nl generate finite volume corrections to the hypercubic symmetry. Only the cubic 
symmetry is exact. We define cubic invariant quantities and compute their average over the cubic 
group. We have thus a set of measures for every orbit of the cubic group, labelled by 



\i=1.3 
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where m = 2, 4, 6. 



A first kind of artefacts that can be systematically cured 12 24 25 are those due to the breaking 



of the rotational symmetry of the Euclidean space-time when using an hypercubic lattice, where 
this symmetry is restricted to the discrete hypercubic H{4) isometry group. However, as already 
mentioned, we have also finite volume effects which break H{4). We therefore need to adapt the 
method. One idea could be to generalise it to a cubic symmetry. This happens not to be practical 
due to too few cubic symmetric orbits for a givenp^. We choose another approach motivated by the 
fact that the lattice action is indeed -ff(4) symmetric and that finite volume effects are expected 
to be small at large momenta compared to finite lattice spacing artefacts. We therefore use a 



slight variation of the method described in 24 ,25 .■ we apply it to the cubic orbits of Eq. (3.11) 



keeping track of p4 which is not an H(A) symmetric quantity. 
Defining the H{A) invariants 



p 



[4] 



p 



4 



pi 



(3.12) 



it happens that every cubic orbit Eq. (3.11) has a well defined set of values for these H(A) invari- 



ants, but several cubic orbits may have the same -ff (4) invariants. We will neglect p'®^ which plays 
no role on small lattices. We can thus define the quantity Zg{ap^) averaged over the cubic orbits 

as 



ap4, a Aqcd 



(3.13) 



We expect the hypercubic effects to be O(a^) lattice artefacts and therefore to be expandable 
into powers of a^. This would of course trivially be the case if a^p^ <^ 1 since then, for example 

However, aiming at 



a2p[4]^p2 <, a'^ p"^ <^ I (xve take on purpose this quantity which will be seen to be dominant). 
Then a Taylor expansion of Eq. (3.13) will ensure the artefact to be O(a^) 

measuring Zg at large momentum we go up to p"^ ~ 3 — 4. We will assume, and then check, that 



the Z^^** in Eq. (3.13) can be Taylor-expanded around p'^^ = up to e significantly larger than 1: 



ap4, o^Aqcd 



+ 



^hyp.corrected | 



a p ,ap4,a Aqcd) 
a^— + ■■■ 

pA 



where 



/?(aV, a^A^co) 



de 



(3.14) 



(3.15) 



Of course terms proportional to p^^\ p'^^^. 



etc. can be added analogously to the formula, as well 
as terms breaking H{A). However we have found that our data were not accurate enough to allow 



fitting them, and that using only Eq. (3.14) and Eq. (3.15) gave satisfactory fits. Now we must 



describe how we fit the functions appearing in the r.h.s of Eq. (3.14). 
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3.3.1 The sliding window fit (SWF) 



We consider all values of a?p^ in a range: a^p^j^ < a^p^ < a^p^nioi im ^ach of which contains a 
set of cubic orbits. We choose an integer width w (we will use w = 10 in numerical applications) 
and define a window as the set of 2w + 1 values of a^p^ around a a?p1^^ter ( contiguous values 
below a^pl^ater ^^"^ ^ many above). There are as many windows as values of a^Pcenter such that 
all of them are in the range [o^Pmin in; "^^Pmax in]- This defines the range of interest o^Pmin out — 

^ Pcentcr — ^ Pmax_ouf 

For every window we use for the fit all cubic orbits corresponding to the values of a^p^ in 



the window. We fit, according to Eq. (3.14), 2w + 2 parameters which are the 2w + 1 values 
of Z^yp-™'''''"'*'"^(a^p^, a^AQ(-,p) within the window, and one common value of /2(a^p^g„(g^, o^Aq^^q). 
The dependence in these parameters is linear and thus the fit amounts to invert a matrix. It is 
clear that for any a?p^ the Z^^P-™"''*^*'^'^(a^p^, o^Is^q-q) is fitted every time a^p^ is within a window, 
i.e. 2 w + 1 times. We keep as the final result only the result of the fit when a^p^ is the center of 



the window. At the end of the fit, for every a^pl^^^^T. in the range a^Pmin out — ^"^Pcenter — ^"^P"^ 



max_out 



we have, as expected, a fitted value for both functions of the r.h.s of Eq. (3.14). 

We can then study the function i?(a^p^, a^AQ(^j3). As will be reported later (see for instance 
Fig. |4]) we find that a reasonable approximation for R is 

R{0^P^, O^A^Cd) = Ca2p4 + Ca4p4 • (3-16) 

This leads to the one window fit. 



3.3.2 The one window fit (OWF) 

We tune w such that only one (or at worst two) windows are included in the range [a^Pmin in; '^^Pmax ii 
We then perform the fit for that window according to the equation 

ylatt/ 2^2 „4^[4] „6^[6] ^2 \2 \ _ yhyp^corrected / 2^2 „2a2 \\„ „2P 

+ Ca4p4 aV"^ • (3.17) 

This fit gives 2w + 3 parameters which are Z^^^-'^°"'^'^^^'^{a'^p'^ , a"^ Aqqj^) for all a^p^ in the window, 
i.e. in the range [a^p^in in' '^^Pmax in]' (^^ the range size is even one value is eliminated) and the 
parameters Ca2p4 and Ca4p4,. 



3.4 Other lattice artefacts 

There are ultraviolet artefacts which are functions of p^ and are thus insensitive to hypercubic 
biases and not corrected by the above-mentioned method. They will be corrected simply by 
assuming a term linear in a^p^ in the final fit and check that the coefficient scales correctly for 
different lattice spacings. 
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To take into account the space-time anisotropy, which is a finite volume artefact, we can check 
the dependence of on a anisotropic quantity such as P4 — We did not see any sizeable 

effect of this parameter. 

Finite volume artefacts are also studied as usual by a comparison of runs at different volume. 
We expect a small effect at large momenta and our checks confirmed it, as well as the analysis 
of (TI. We will not consider this artefact anymore. 



4 Lattice results and hypercubic corrections 

The hypercubic artefacts generate on the raw lattice data, Z^^^^ ^ the so called "half-fishbone" 
structure 12 shown in fig[lj In this figure all the points labelled as explained in Eq. (3.10) and 
Eq. (3.11) are plotted. The color code shows the value of the ratio which is between 



0.25 and 1. The closer to 1 means the less "democratic" or "tyrannic" ones. We see, as expected, 
that the tyrannic points are more affected by the artefacts. We also see that the gap between 
_^iatt^^2p2^ ^4p[4]^ ^6p[6]^ ^p^^ ^ given ij? cau be as large as 0.07, i.e; about 10%. Taking 

a naive average, without a correct treatment of this artefact would leave a systematic upward 



shift of about 5 %. In section 3.3 we have developed a non-perturbative method to cure this 
artefact (this method being indeed splited in two options, the SWF and the OWE). There exist 
two other methods. The oldest one is the "democratic selection". It amounts to keep only, say, 
the points with ratio < 0.3 in fig [Tj One sees that it eliminates a lot of data and it is still far 
from the "egalitarian" result (the lowest curve in fig [T]) which we consider as better. We will thus 
not consider further this democratic selection. The second other method to correct hypercubic 
artefacts uses perturbative calculation and is detailed in the next section. 



4.1 Perturbative correction 



The perturbative method 26 consists in computing at one loop in lattice perturbation theory 27 



then, assuming that the lattice spacing artefacts are reliably described by the 0{g'^a'^) terms 
thus obtained, in subtracting them from the lattice data. This method has been applied to 
quark bilinear operators in [2^. Eor comparison we have applied here this method following the 
prescription described in section 3.2.2 of 28 and also a variant of it. 



Eq (24) in [26] may be written in Landau gauge as 
Zf^*(aV) = 

Cg3 



where, to follow [28], = gl^^^^^J{127r^), with gl^^,^^^ 



oV) + 9^{bqi + Cgi aV + Cq2 log(aV) + 

(4.1) 

/< plaquette >. The coefficients 



fl2 [4] 

+ ^94^:^log(a p )) , 
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.9 =< ratio i — i — 

.8 =< ratio < .9 ^ 

.7 =< ratio < .8 
.6 -< ratio < .7 c 
.5 -< ratio < .6 
.4 -< ratio < .5 

.3 -< ratio < .4 :-- •- 
.2 =< ratio < .3 - - 

, hyp corrected OWF ; * 



Figure 1: This plot shows the raw data for /3 = 3.9, ^^^"(a^p^ a^pW, a^pt^l, ap4, a^A^j^.^) in Eq. (3.14), 
in terms of a'^p^ in the horizontal axis. The "half-fishbone structure" due to hypercubic artefacts is 
clearly seen. There is one point for every cubic (3-D) orbit. The color code classifies the data according 
to their degree of "democracy" measured by ratio = The lowest plot corresponds to the 

non-perturbatively hypercubic corrected data (or "egalitarian result" ) resulting form the one-window fit 



hyp_corrected / 2 2 



features appear for all /3's. 



(a^p^, o^Aq^^j-)) in Eq. (3.17). The data correspond to (3 = 3.9, o/i = 0.004, but the same 



are defined using the notations of eq. (24) in 26 



.(2,4) 



.(2,1) 



80 10 



^"^^ 240 2 60 ' 



101 11 

120 ~ 30 



Co . 



4.1.1 Prescription with 



Using the prescription of eq. (35) in 28 , for every cubic orbit we define the substracted quantity 



as : 



~2 f 2 ~2 2 '^1 1 / 2 ~2 \ Ct p Qj p ^ 

g \ Cgi a p +Cg2a p log(a p ) + Cg3—^ + Cqi—^ log( 



2 

a p ^ 



(4.2) 



The result of this substraction is plotted in the l.h.s of fig. [2j The half-fishbone structure is still 
clearly visible for a?p^ > 1.6. We then try another prescription. 
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Figure 2: L.h.s: data of fig.[T]corrected by the perturbative subtraction, formula of Eq. (4.2 ), section 4.1.1 



R.h.s: same exercice with the prescription Eq. (4.3) with Eq. (4.2). The color code is the same as in fig[T] 
and ratio = p^^^ / {p^)'^ . The horizontal axis is a^p^ for both. 



4.1.2 Prescription with 



The trace of Eq. (3.9), which introduces p^, had to be applied in eq. (24) of [26] to obtain 
Eq. (4.1 ). We will now expand in before performing the trace and then keep the 0{g'^a'^) terms 



for substraction. This gives, using again 26 



/ fl2 [4] 2 [4] 

~2 / 22 2i/2 2\ f Uj jj ^ y 1 / 2 2\ 

- 9 [Cqi a p + Cg2 a plog(a p ) + Cg3—^ + Cg4-^ log(a p ) 



(4.3) 



where 

<3 = V + ^^''^ - ^ - ^ = + ^ . (4.4) 



-(0,1) Q r c(o,i) 

^''^ 6 80 10 ' 6 



This result is plotted in the r.h.s. of fig. |2} With this variant, the half-fishbone is significantly 
reduced but the dominant 0(4) artefact is overcorrected : a linear behaviour in a^p^ is clearly 
visible, larger than with the first prescription. 

4.1.3 Lessons about the perturbative method 

We see that the two prescriptions start differing significantly at a^p"^ ~ 1, which is not surprising 
since higher order terms become significant, for example, an a*p^*^ is also of order 1 for tyrannic 
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points, while a^p^ — a^p^ ~ 0.3. The perturbative method goes in the right direction, but it is 
impossible to know a priori its quality without performing the tests we propose here. The method 
contains several ambiguities: what to take for the coupling constant ? use or p^l Contrarily to 
the non-perturbative method it provides both the hypercubic corrections and the 0{o?p^) ones. 
Conceptually the perturbative method is very useful as it exhibits qualitative aspects which may 
guide the use of the non-perturbative one, for example it justifies the smoothness of the variation 
of the derivative R Eq. (3.15) as a function of aP'p^ as well as that of the slope in o^p^. Finally we 



shall see that it gives results similar to the non-perturbative ones. 

4.2 Non perturbative hypercubic correction 
4.2.1 Sliding window fit versus one window fit 



We now apply the non-perturbative correction. We only expand in since the higher order terms 



turn out to be negligible in our momentum range. In section 3.3 we have presented two types of 



fits, similar in spirit: the sliding window fit (SWF) described in section 3.3.1 which amounts to 



using Eq. (3.14) combined with Eq. (3.15) and the one window fit (OWE) described in section 3.3.2 
which amounts to using Eq. (3.16). In the l.h.s of Fig. [31 we show, in the case of /3 = 3.9, the 



comparison of hypercubic corrected data after applying OWF and SWF. The difference does not 
appear to be large which is rather encouraging. The OWF gives a slightly smoother result. For 
this value of /3 the chi-squared is not good (see table [2]) but remember it uses only two hypercubic 
parameters. Chi-squared for the other /3's are better. 



4.2.2 Half-fishbone reduction test 

We need also to apply the half-fishbone reduction test as in the perturbative case, i.e. to subtract 
to the raw data of every cubic orbit the hypercubic correction. We present the OWF result. From 



Eq. (3.17) the subtraction amounts to 



Ca2p4 ' 



p2 



(4.5) 



The result is shown in the r.h.s of Fig. |3} one point per cubic orbit. The non-perturbatively 
corrected Z^^^-^°"'^'^^^'^{a?p'^ ^a^K^Q^^) of Eq. (3.17) is represented by the black line in the r.h.s of 
Fig. |3| It is well in the middle of the subtracted points, as we would expect. 

It is seen that the half-fishbones have been strongly reduced. One sees a remainder of these 
artefacts due to the less democratic, or "tyrannic" points. These have only one non-vanishing 
component or one large and a very small one. These points are not so many as seen in the plot, 
their orbits are small which explains the larger error bars. 

We have checked that these tyrannic points have indeed a small impact on the hypercubic 
corrected result. 
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Figure 3: On the l.h.s we compare for (3 = 3.9, = 0.004 the non-perturbatively corrected Zq from the 



shding- window fit (SWF), section 3.3.1 and from the one- window one (OWF), section 3.3.2 To show 
better the running we have also subtracted the 0{o?p'^) artefact which will be computed below. On the 



r.h.s we show, using the OWF, the non-perturbatively subtracted data Eq. (4.5). There is one point for 



every cubic (3-D) orbit. The color code and the definition of the parameter ratio are the same as in fig. [T| 
The black line corresponds to the OWF non-perturbatively corrected result: z]^^^~^°"'^'^^'^'^ {a^p^ , a? A^qj^) 
of Eq. (3.17). The 0{a?'p'^) has also been subtracted. The horizontal axis is a?p'^. 



4.2.3 The slopes in p 



[4] 



The sliding window fit solves for every window a value for the slope R{o?'p'^) Eq. (3.15), i.e. the 
derivative dZ^^^^ /d{a?'p^'^yp'^). This allows for a study of the shape of this function R. In Fig. |4] 
we plot this slope R defined in Eq. (3.15) for the three values of {3. We also plot the equivalent 



slope using the perturbative formula with the p^ prescription for (3 = 3.9, section 4.1.2 R 



l-pert 



Cg3+Cg4 log(a p ), defined in Eq. (4.4) and Cg4 in Eq. (4.2). We see that this perturbative slope is 
in fair agreement with the non-perturbative one, explaining the good elimination of half-fishbones 
in the r.h.s. in Fig. [2j 

The three non-perturbative data in Fig. |4]give the impression to be affine (a constant minus a 



linear term) over a rather large momentum interval. This is what is expressed in Eq. (3.16) from 
where we have deduced the one-window fit: a fit over the full range [0.5-3.5] with two hypercubic 
parameters only |^ 

The fitted values for Ca2p4 and Ca4p4, from the one window fit are given in table |2] as well as 
the same divided by g^, since perturbation theory expects at least for Ca2p4 to be oc g'^. Before 



*0f course there are additionally as many hypercubic insensitive parameters as there are values of in the 
range, which are simply the values of Z^''P-™™'^*'"^(a^p^, a^AQ,^^). 
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Figure 4: The three lattice slopes R defined in Eq. (3.15) for the three lattice spacings and the pertur- 
bative one Rpert = + Cq4 log(a^p^) for (3 = 3.9, in terms of a^p'^ in the horizontal axis. 



(3 




Ca2p4 


Ca4p4 


Ca2p4/ 9^ 


Ca4p4/ 9^ 


xVd.o.f 


3.9 


0.00689 


0.067(4) 


-0.0149(10) 


0.044(3) 


-0.0097(7) 


4.1 


4.05 


0.00456 


0.065(3) 


-0.0144(5) 


0.044(2) 


-0.0097(3) 


0.53 


4.2 


0.00303 


0.055(11) 


-0.0124(4) 


0.039(8) 


-0.0089(3) 


0.98 



Table 2: Results for the slope in a^pW/p^ and a^pW and the same divided by g"^ in the one window fit. 



dividing by 9"^ a small scaling violation is apparent, which corresponds to the non overlap of the 
curves at different /3's in Fig. 
the scaling. The in table |2 
lower end of the plot. 



4| It appears on table [2] that dividing by 9^ improves significantly 
is not good for (3 = 3.9, apparently due to some structure at the 



5 Running including (A^) corrections from OPE 

In this section we will check the running of Zq. To this purpose, we shall use both the formula 
Eq. (A. 24) derived in the appendix [A| to which we add a lattice artefact term oc a^p^ not yet 
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gH2 fron one ulndou tree ' — 1 — i 
average of gfi2 fron OHF 
gHS fron slidinE ulndous tree ' — *t — i 
average of gfi2 fron SHF 
£02 fron one ulndou CM - 
average of OHF uith CM 


1 


gH^ tron s±j.aj.ng uj.naous lm 

average of SHF uith CM 

I'ri2 fron the fit on nerged data 
gHS f ron the strong coupling constant 

t 























ca2p2 fron one ulndou fit (OHF) > 
average of ca2p2 fron OHF • 
ca2p2 fron sliding ulndous fit (SHF) ^ 
average of ca2p2 fron SHF ■ 



Figure 5: We plot the values of the fitted slopes Ca2p2 (r.h.s) and the condensate (A^) (l.h.s) as 
extracted from the l/p^ contribution to the fit, see table [sj |4] and [?[ In the left plot we show the results 
from the OWF and from the SWF. It can be seen that the 0(c?) formula for the Wilson coefficient 
computed by Chetyrkin and Maier [9j of (^4^) (indicated by the "CM" initials) is about 20 % below 
the tree level result. We show the value obtained from the merged results of the three lattice spacings, 
table [7j Finally, for the sake of comparison we show the result from the strong coupling constant of [?]. 
The horizontal axis is in fm^. 



subtracted: 



1 + 



+ Ca2p2 0?P^ ■ (5.1) 



and a formula including only an OPE correction with a tree-level Wilson coefficient, 

yhyp.corrected/ 2 2\ r^pert RI' / /2\ RI' ( ( 1 , Cio^erp2 ^ 2 2 

l« P J = l/^ )Cqz^ \y2'^y^^) j ( ^ ~2 ) + Ca2p2 « p 



(5.2) 



where Cioverp2 = 5'^(^^)/12. We use /x' = = 10 GeV as the renormalisation scale ; c^^^{p'^ , iJ?) is 
computed from the four loop perturbative running of Zg [l3| : 

^pertRI'(-p2^^2^^J 

CoZg(P ) = „pcrtRI'/ 2 2 N ' ^-^'3) 

[f^ 1%. 



are/ 
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ZQ rron one Mindou fit <OHF> 
linear fit of ZB fron OHF ■ 
ze fron sliding Mindous fit <SHF> 
linear fit of ZB fron SHF ■ 




Figure 6: We show the value of Zq"^^ defined in Eq. (5.1) for all three lattice spacings as a function of 
the bare coupling constant g'^ = 6//3. 



c^^(p^,/i^) is the three loop Wilson coefficient of {A^) in the expansion of Zg |9j and the ratio 

^(/^)) 1 - 0.1317 _ 0.5155 



9 V M 



1 - 0.1317 a^{p) - 0.5155 a^{p) 



(5.4) 



was obtained in appendix A. We express the lattice spacing (cut-off) dependence as a dependence 
in (7bare- Z^'^''* (/^^! fi'bare) the pcrturbative contribution to Zg at the scale /i in the RI'-MOM 
scheme. In other words, 



1 + 



arcy 



.RI' _ 



a{^^)] c^l g,a(/i) 



32 p2 



(5.5) 



From now on ^p"''* will refer to Z'p®''* We fit three parameters: Z'^'^'^^, Ca2p2 and alternatively 
Cioverp2 (which amouuts to a tree level treatment of the C2Zq coefficient) or the vev g'^{A'^). In order 
to estimate the systematic errors we will treat in parallel the one window fit and the sliding window 
one. The results are reported in table |3| table |4], Fig. [5] and Fig. |6} The coefficient Ca2p2 obviously 
refers to an 0(4) invariant lattice spacing artefact which is not detected by our non-perturbative 



20 



/3 




^pcrt 




9^ ^A^) tree 




3.9 


0.00689 


0.726(5) 


0.0201(13) 


3.20(38) 


2.62(31) 


4.05 


0.00456 


0.742(5) 


0.0200(15) 


3.09(65) 


2.57(54) 


4.2 


0.00303 


0.760(3) 


0.0194(8) 


3.23(55) 


2.74(47) 


average 






0.0201(3) 


3.18(28) 


2.64(23) 



Table 3: Results for '*(10GeV) and Ca2p2 Eq. ( |5.l[ ) from the one- window-fit and the estimated g^{J^) 

' ' 7pert 



vev from the 1/p term and from the the Chetyrkin-Maier [9| (CM) Wilson coefficient. Notice that Zq 
and Ca2p2 from these two fits are very close. 



/3 




^pert 


Ca2p2 


g^{A^)tree 


g\A^)cu 


3.9 


0.00689 


0.741(3) 


0.0161(9) 


2.07(37) 


1.70(31) 


4.05 


0.00456 


0.753(5) 


0.0168(14) 


2.13(52) 


1.78(43) 


4.2 


0.00303 


0.771(3) 


0.0164(9) 


1.59(60) 


1.36(51) 


average 






0.0165(6) 


1.99(26)(27) 


1.65(22)(27) 



Table 4: The same as in table |3] using the data from the sliding- windows-fit to hypercubic corrections. 



hypercubic correction method. We see in the right plot of fig. |5] as well as in the tables that this 
coefficient scales very well when expressed in lattice units, as it should be. The coefficient of 1/p^, 
if it is related to a vev (A^), should rather scale in physical units. We see in the left plot of fig. [s] 
that a constant value is rather well verified although with large errors. The results presented in 
table [3] and [4] show that the estimates for g'^{A'^) from OPE expressions with Wilson coefficient 
from the Chetyrkin-Maier (CM) three-loop expression are systematically about 20% below the 
ones from the tree level one. 



5.1 Analysis from the non-perturbative hypercubic corrections 
5.1.1 Comparison of the running from the OWF and the SWF 

From tables [3] and [4] we see that {A^) is systematically larger for the OWF than for the SWF. At 
first sight it seems surprising since the OWF and SWF hypercubic corrected data are very similar, 
see Fig. [3J One reason is the correlation between Ca2p2 and {A'^): Ca2p2 is also systematically 
larger for the OWF than for the SWF. This correlation is understandable as the a^p^ increases 
with p"^ while l/p^ decreases. This is compensated by a Z^"^^^ smaller for the OWF than for the 
SWF. We will consider these differences as a systematic uncertainty in our fits and count them in 
the errors. 



5.1.2 Dependence on the fitting range 

An additional test is to look for the effect of the fitting range. The results are shown in table [5] 
and Fig. u\ One sees again a correlation between the Ca2p2 slope and g"^ {A"^). Both decrease when 
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Figure 7: In the r.h.s plot we show how g'^{A'^), fitted with the CM Wilson coefficient, depends on 
the upper range of our fits, always starting at a^p^ = 0.5 for /3 = 3.9. The points correspond to 
a^p"^ < 2.0,2.5,3.0,3.5. The r.h.s plot also shows the a^p^ artefact slope. We find again a positive 
correlation between both series of data. The l.h.s of the plot illustrates this, showing for the same data 
how the fitting function depends on the upper bound of the range. 



the fitting range shortens while correlatively Z^'^'^ increases. The l.h.s of Fig. [7] explains how this 
happens: when the range is shorter, the error bars allow for a less curved fit. But never the fit 
reaches a value such that (A^) disappears. The shortest window 0.5 < a'^p^ < 2.0 gives the 
smallest value for g'^ (A^) but still 4 sigmas away from 0. 



5.2 Analysis from the perturbative hypercubic corrections 

It is then useful to check if similar results are obtained after a perturbative correction to the 
hypercubic artefacts has been applied. We have used two prescriptions to apply the perturbative 



corrections. With the data obtained from the prescription, section 4.1.1 , we perform an average 
on all the cubic orbits of every p"^ after a democratic selection < 0.3. This leaves us 

with not too many points and it results in rather large statistical errors. We then perform the 
same running fit than on the non-perturbatively hypercubic corrected results: we fit with one 
perturbative running contribution, one l/p"^ contribution and one oc a^p^ artefact. With the data 



from the prescription, section 4.1.2, we perform the same fit over an average on all the cubic 
orbits of every p"^ without any democratic selection, since the hypercubic artefacts have already 
been efficiently reduced. The results are in table [6j 

A first remark is that the Ca2p2 coefficients are compatible with zero, indicating that the 
perturbative correction has efficiently eliminated this artefact. The coefficient of the 1/p^ non- 
perturbative contribution is found different from zero, in the same ballpark as the results from 
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Upper bound 


^pcrt 


Ca2p2 






2.0 


0.754(6) 


0.0089(21) 


1.58(39) 


1.28(32) 


2.5 


0.745(6) 


0.0130(18) 


2.05(37) 


1.67(30) 


3.0 


0.733(5) 


0.0175(15) 


2.73(36) 


2.22(30) 


3.5 


0.726(5) 


0.0201(13) 


3.20(38) 


2.62(31) 



Table 5: /3 = 3.9: results for the Zf" (^= 10 GeV) and Ca2p2 Eq. (pdj) from the one-window hypercubic 
corrected data and the estimated g^{A^) vev from the plotted as a function of the upper bound of 

the fitting range (in GeV). 



prescription 


^pcrt 


Ca2p2 


g'^{A^) tree 






0.712(11) 


-0.0026(23) 


2.98(1.49) 


2.53(1.23) 




0.745(3) 


-0.0061(8) 


1.76(35) 


1.45(29) 



Table 6: For /3 = 3.9, results for the Zl'''\fi'^ , gl^^^) {fi= 10 GeV) and Ca2p2 Eq. ([5^ and g'^{A'^) from 
the lattice data after a perturbative hypercubic correction. g^{A^) is estimated at tree level from the 
contribution. 



the non-perturbative hypercubic correction, in tables |3] and |4j The prescription has too large 
errors to be conclusive but the is five sigmas away from zero, very similar to the results in 
table 4j The value of Z^'^^^ for the prescription is rather low but compatible within less than 
two sigmas from the result for {3 = 3.9 in table |4] 



5.3 Running of 

It is interesting to consider the dependence of Z^^'^^ as a function of g'^. This is plotted in fig. [6] 
both for the OWF and the SWF. It is strikingly linear, specially for the OWE. Indeed from [26] 7 
eq. (24) perturbation theory gives a linear dependence with a slope ~ —0.19. This comes from 



the coefficient bqi in Eq. (4.1) which is multiplied by g^. In our case we find from the OWE 

-- 0.737(3) - 0.313(6) (^Le - 1-5) , 
- 0.766(3) - 0.324(6) (^Le - 1-5) ; 



Zf*((10GeV) 



Zf'^*((2GeV)^ 



2 ^2 N 

' i/bare/ 
2 \ 
' dhare) 



and from the SWE 



Zp((10GeV)^^7Le) 
Zr\{2GeVr,glJ 

We see that the coefficients of g^, 

dZ^^^\{2GeVf,gl 



0.751(2)(7)- 0.273(6) {^Z%s) (gL. 



0.780(3)(7) - 0.284(6) ( 



0.002 \ ( 2 
O.OiO) li/bare 



1.5) , 
1.5). 



dg' 



-0.329(6) OWE 
-0.287(26) SWE 



(5.6) 



(5.7) 



(5.^ 
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are significantly larger than the perturbatively expected, -0.19. But the linear behaviour predicted 
by perturbation theory is well verified, especially for OWF. 




Figure 8: The merged plot with the OWF results at /? = 4.05 and /3 = 4.2 rescaled to the /3 = 3.9 thanks 
to the ratios of Zq^^^ given in table [sj The l.h.s shows the data corrected for all lattice artefacts. The 
r.h.s shows the same data furthermore corrected by the perturbative running factor up to 10 GeV. The 
horizontal axix is in GeV ^. The black line on the l.h.s corresponds to the global fit with perturbative 
running and CM (three loops) wilson coefficient for the l/p"^ term. The black line on the r.h.s corresponds 
only to the times the three loops wilson coefficient added to Zq'''^*((10 GeV)^, 6/3.9) = 0.726 



5.4 Merging the three lattice spacings. 



From Eq. (5.1) and Eq. (5.5) it is clear that 

rvKV (2 2 \ ^hyp.corrected/ 2 2\ ^ ^2 2 /r q\ 

In this section we use the one window fit, section 3.3.2 and the momentum is now expressed 



in physical units. For the coefficient Ca2p2 "we use the values in table [3] Ca2p2 = 0.0174. The three 
^Ri ^Qj, ^]^g three /S's do not match due to the running of Zg as a function of the lattice spacing. 
To make them match it turns out that it is enough to take into account the ratios of Z^^^^^s given 
in table |3| We plot on the l.h.s of fig. |8]the three sets of data where the /3 = 4.05,4.2 ones have 
been rescaled to the P = 3.9 scale. We see a rather good overlap. There is however a flattening at 
the right end of every /3 which stays within one sigma from the other /3's. We understand it as a 
failure of the hypercubic artefacts treatment. On the r.h.s side of fig. |8]we plot the same number 
corrected for perturbative running by the multiplicative factor 0.726/ZP'^''*(p^, 6/3.9), where 0.726 
is taken from table [7| The black line is just the non-perturbative contribution added to 0.726. 
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The comparison of both plots in Fig. [8] is enhghtening. We see that the non perturbative term 
contributes about one half of the change between the smallest momenta and the largest ones. 
Both the perturbative running and the non perturbative contribution are convex, which makes it 
difficult to disentangle them. But we also see that the perturbative running cannot account for 
the full variation of the data. The best-fit parameters resulting from this merged analysis can be 
found in table [7| where we use g^.^^^ = 6.0/3.9 since we have rescaled all the data to the /3 = 3.9 
one. The values in table [t] for g'^ (ii^) {A^) ^2 merged turn out to be rather central in the set of values 
of tables and El The value for ZP'^''*(10 Gel^, 6./3.9) is in very good agreement with Eq. (5.6): 
ZP^'"*(10Ge 1^,6.73.9) = 0.725(3). 



^pcrt 


9'^ {A)tree 




0.726(2) 


3.13(43) 


2.55(36) 



Table 7: Merged data from three /3's: results for the Zq"^^^ {fi= 10 GeV) rescaled to /3 = 3.9, from the 



one- window hypercubic corrected data (OWF) with tree level and the three loop formula, Eq. (A.24). 



5.5 Summarizing 

Many of our results for g'^{A^) are shown in the l.h.s of Fig. [s] We did not plot the range 
dependent data to not overload the plot but they faU in the range covered by the data plotted 
in Fig. |5] To summarise our results and estimate the systematic uncertainty we consider the set 
of values in tables |3, 4][5 and[7j We will make a separate average for the tree level data and the 



0{a'^) (CM) ones since the comparison with estimates from other quantities, such as the coupling 
constant, need to be performed in the same scheme, expansion, order and scale. The scheme is 
MS and the precise implementation is detailed in the appendix. We computed an average of all 
above-listed data, weighted by their inverse squared error. The inverse squared statstical error is 
the sum of the inverse squared errors. The systematic error is taken such as to incorporate the 



central values within the error-bars. This is rather conservative. For Z^'^''* we average Eq. (|5.6|) 



and Eq. (5.7) with a similar method. We get: 



g\fj?){A^),2,,,, = 2.45(14) GeV^ /i = 10 GeV , 

g\l?){A^),.CM = 2.01(11) (1°:?^) GeV^ /i = 10 GeV , 

Zp((10GeV)^^?Le) = 0.744(2)(7)- 0.311(6) (^^Lc - 1-5) , 

ZP-H(2GeV)^^7Le) = 0.773(3)(7) - 0.323(6) {t'oZ) idLe ' 1-5) ; (5.10) 

where the first error is statistical and the second the systematic one. The values of Z?-^' {p'^ , g'^^^.^) 



may then be derived from Eq. (5.5): 



Zf{{2GeVr,glJ = 0.805(14) - 0.336(6) (;n42) (^^Le " 1-5) • [CM] (5.11) 
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The results obtained from the perturbative hypercubic correction, table |6] have not been used in 
the final estimate (because we did not compute them for all P's) but they fall within the bounds 
at less than one sigma. 

Finally two lines in the l.h.s. plot of fig. [sjshow the results for g'^{A^) obtained from those in 
tab. 3 of ref. [T]: 



9' {A') 



10 GcV 



4.1 ± 1.5 GeV' 
2.5 ± 0.9 GeV^ 



leading log 



(5.12) 



These values in Eq. (5.12) come from a totally different quantity: they have been extracted from 



the running of the ghost-gluon couphng constant. The results of ref. |7j were obtained by applying 
an OPE formula including a Wilson coefficient approximated at the leading logarithm and at the 
order (9(a^), but expanded in terms of ax- Then, in order to be properly compared with the 



results of this work, the values of Eq. (5.12) incorporate the correction by the effect of expanding 
the OPE formula in terms of the running coupling in MS. The lattice spacing applied in ref. ^ to 
get a physical scale, a(3.9) = 0.0801 fm, was also slightly smaller than the one used in this work 



(see, for instance, tab. [T]) and this has been also taken into account in obtaining Eq. (5.12). 

As a matter of fact, Eq. (5.12) exhibits a slower convergence of the perturbative series of 



the Wilson-coefficient as in the present paper. From tables [sj |4| |5] and [T] we see that the O(a^) 
estimate is about 20 % below the tree level while in table 3 of |7j it is about 45 % below the 
leading logarithm one. The two estimates agree rather well within the present accuracy. 



5.6 Conversion to MS 

The conversion to MS can also be performed from Z^^"^^ or Z^^' which contains a non-perturbative 
contribution. Usually, in the literature, the values are assumed to be perturbative. The conversion 



of ZP'^''* into MS will use the standard perturbative conversion formulae 13 . We get 



^MSpert((2GeV)^^Lc)/^^((2GeV)^^Le) = 0.97 , 
P-*((2GeV)^^Le) = 0.750(3)(7) - 0.313(20) (^Le - 1-5) • (5.13) 



The central value of this result is about a 2 % systematically below the results of 28 . This can 
be presumably interpreted as a small systematic correction due to our subtraction of the non- 
perturbative contribution in obtaining Z^*^"^^ in the RI'-MOM scheme and converting it to ^^Spert 
at 2 GeV. 



From Eq. (A. 24) it is easy to see how to include the g'^{A^) non perturbative contribution. 
Up to now we have applied the result of the appendix using RI'-MOM for the perturbative part 
and MS for the ratio in the corrective parenthesis. Had we wished to use the MS scheme for 
the perturbative contribution, we would have the main inconvenience of not to know the OPE 
contribution for Zq defined in the MS scheme. However, only in the aim of roughly estimating the 
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non-perturbative correction, we can assume the OPE corrective parenthesis to remain the same 
and then get 



^MSnon-perturbative^^2GeV)^^?Le) = 0.781(6)(21) - 0.326(21) (^Le ' 1-5) 



(5.14) 



Notice that the non perturbative contribution is about 4% at 2 GeV. Nevertheless, the results 
of [28] have been obtained at momenta larger than 2 GeV. Although their estimates and our 
Z^^nay agree with each other, a subtraction of the non-perturbative contribution in obtaining 
^pert^ probably still required at the fitting window of [28] but not applied, could explain the small 
discrepancy of about a 2 % that we discussed above. The discrepancy is anyhow only affecting 
the conversion of Z^'^' to Z^p^=^* at 2 GeV. In Table 6 of 
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one finds results for obtained 



at 2 GeV by applying the standard perturbative conversionformulae 13 . Indeed, they turn out 



to be in between our results in Eq. (5.13) and Eq. (5.14). Quoting in order Eq. (5.13), ref. 28 



and Eq. (5.14) we get: 0.738, 0.751 and 0.769 for /3 = 3.9 and 0.755, 0.780 and 0.786 for (3 = 4.05. 



5.7 Comparison of different estimates for the gluon condensate 

Let us now compare the present estimate of g'^{A^) to previous ones at Nf = 2 and Nf = 0, all 
taken at the renormalisation scale of 10 GeV and, when needed, transformed to the very precise 
renormalization scheme for the OPE expansion defined in the appendix |Aj 

In j4] a quenched study of Zg using Wilson-Clover and overlap fermions ended with values of 
(A^)mom in the range 2.67-3.2 GeV^ with typical errors of 0.3 GeV^. Notice that this computation 
was performed only up to leading logarithm for the Wilson coefficient and that the choice was to 
expand the perturbative series in terms of the coupling renormalized in the MOM scheme (this 
is why we use for the VEV the label MOM). Then, we can apply the expressions derived in the 
appendix [a] to obtain the estimates for g'^{A^) in the above-mentioned renormalization scheme 
appearing in table [s} However it is advocated in [i] that the l/p^ contribution only increases by 
10% when going from MOM to MS. On the other hand we have seen that between tree level and 
three loops a decrease of 20% was observed. In [i] an artefact oc a/p^ was observed. We do not see 
it in the present analysis since the scaling of g'^{A^) as a function of the lattice spacing indicates 
no visible contribution dependent on a. 

In [6] a summary was performed of different estimates of g'^{A^) from gluonic quantities at 
Nf = : from the three gluon vertex with equal momenta on the three legs (symmetric) 
and from the three gluon vertex with one vanishing momentum (asymmetric), the ratio between 
the ghost and gluon propagators and from the ghost and gluon propagators, using Taylor's 
theorem. The ones involving gluon and ghost propagators agree fairly well but the latter is the 
most accurate. It gives gq-^A^) = 5.1^°;^, although the applied OPE formula was obtained by 
expanding the involved perturbative series in terms of ut- After the appropriate transformation, 
one obtains the result shown in table [s} We also quote in the table the estimate of g'^{A^) from 
the symmetric three gluon vertex, more precise than the one coming from the asymmetric vertex, 
and that appeared to be much higher than the estimate from ar and compatible with that from 
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quark propagtor. In the case of the three-gluon estimates, no available O(a^) Wilson coefficient 
can help us to go beyond the leading logarithm approximation. However, either comparing the 
leading-logarithm estimates of the ones approximated at the order O(a^), a clear discrepancy 
(by a factor of about two) appears between the estimates from ghost and gluon propagators and 
those from vertices or the quark propagator. This could imply that some systematic uncertainty 
is not completely under control. One might, for instance, guess that 1/p^-contributions can be 
invoked to reduce that discrepancy. For this to happen, the 1/p^-contributions had to be negative, 
and to tend to increase the estimate of g'^{A'^), for the OPE formula of a^; while it had to be 
positive, and reduce g^{A'^), for the quark propagator. Indeed, although no stable fit including 
l/p'^-contributions can be performed, the sign seems to be the right one for ax in |7|. Also the 
right sign of the contributions to Zg is found in ref. [il. 



measurement (GeV^ 



Nf 


order g^{A^) 


z. 


ar 


3-gluon 





LL 


9.4(3) 


5.2(L1) 


10(3) 




9.0(3) 


3.7(8) 




2 


LL 


2.7(4) 


4.1(L5) 






2.55(36) 


2.5(9) 





Table 8: Comparison of estimates of g^{A^) from different quantities an Nf = and Nf = 2. All are 
taken at the scale = 10 GeV. LL means leading logarithm for the Wilson coefficient. O(q^) refers to 
Chetyrkine and Maier computation. 

In [T] the strong coupling constant was computed along similar lines to what is done here, on 
the same set of ETMC gauge configurations with Nf = 2. The necessity of a non-perturbative 
oc l/p^ contribution was also found and the resulting condensate, 5'^(^^)ioGeV = 2.3(8), obtained 
through an OPE formula approximated at the 0(a^)-order and expanded in ax, can be properly 
transformed to give the value of table [sj also quoted in Eq. (5.12), which agrees strikingly well 
with the result of tab. [T] (the one we also quote in tab. [8| or that of Eq. (5.10 ). The value obtained 

In [Tol Martinelli 



through a leading-logarithm-approximated formula is also displayed in tab 
and Sachrajda proposed a criterium to validate the use of operator expansion which we apply in 
this paper. They concluded that one should compare the difference of the highest order of the 
perturbative expansion for two different quantities with the non-perturbative contribution and 
check that the former is small compared to the latter. We have compared the highest order of the 
perturbative expansion of Zg with the l/p^ contribution and find that the ratio ranges between 
1/10 and 1/3 depending on the momentum. This is a good indication of the validity of our use of 
the operator expansion. Had we used the perturbative expansion only up to 0{a) the criterium 
would not have been fulfilled. 

Furthermore, all these estimates in table [s] show a clear tendency of a decrease of g^{A^) from 
Nf = to Nf = 2. This might support an interpretation of g'^{A^) as originating in instantons 
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'We have taken into account the different lattice spacing in j^. 
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since the instanton density should decrease with hght dynamical masses. 



6 Conclusion 

We have studied with care the twisted quark propagators produced on the ETMC set of Nf = 2 
gauge configurations. Our goal was to concentrate on two major issues: the correction for lattice 
spacing artefacts, particularly the hypercubic ones, and the presence of a sizeable non perturbative 
contribution of the operator. The latter is expected to be sizeable because it was seen in the 
quenched case |4|, it was seen in the unquenched study of the strong coupling constant [7j and 
since the Wilson coefficient of g'^{A^) is not small in Zq [Oj. This is an important issue since from 
our estimates it gives a ~ 4% contribution at 2 GeV. A reliable estimate of this non-perturbative 
correction needs a large enough fitting range, which allows to distinguish a contribution from 
perturbative logarithms. But the fitting window is restricted below by infrared effects and above 
by lattice spacing artefacts. We thus need to improve our control on dominant lattice spacing 
artefacts which are of two types: hypercubic ones and oc a?p^ ones. 

Concerning the hypercubic artefacts, we have summarized the non-perturbative correcting 



method 24,25] which we compared systematically with the perturbative results of 26 . Zq has 
very large hypercubic artefacts which display, as a function of a "half-fishbone" very far from 
a smooth curve (see Fig. [T|. We check carefully how these fishbones are "swallowed" by the 
corrective methods. It is worthwhile to emphasize that the "democratic" method, prescribing for 
instance a cut on p^^^p^ to drastically reduce the number of allowed hypercubic orbits, is not good 
enough to eliminate the fishbones and to leave us with a smooth curve for Zq. 

The perturbative method to correct hypercubic artefacts suffers from some options left: what 
to take for the coupling constant, use of or = sin(ap^) ? We tried first to stick to the 



prescription of 26 and use the boosted couphng constant. This reduces the hypercubic artefacts 
only up to a^p^ ~ 1.6 (see Fig. [2| l.h.s). Guided by the test on the fishbone reduction we then 
propose a prescription based on the same perturbative formulae but using p^. For Zq this reduces 
the hypercubic artefacts up to a^p^ ^3.5 which has been our upper bound in this work (see Fig. 
|2) r.h.s). 

We test also the non-perturbative method to correct hypercubic artefacts. We use two pre- 
scriptions. The first one uses a sliding window and the second one uses only one fitting window 
on the full momentum range. 

We find that the hypercubic artefacts are sufficiently well described and cured by two terms: 
oc a^pW/p^ and oc a^pt^'. We fit the coefficients of these quantities and check their scaling with (3. 

From the resulting hypercubic corrected function Zg(a^p^, a^Aq^j-,) we perform fits which in- 
corporate the perturbative running, a non perturbative term (presumably related to g'^{A^)), 
and a hypercubic insensitive lattice spacing artefact proportional to a^p^. The fits are good. The 
a?p^ term scales almost perfectly in lattice units, as expected. The g'^{A^) term scales rather well 
in physical units as expected. The accuracy on g'^{A^) is reduced by some correlations in the fits: 
we see a correlation between the method used to correct hypercubic artefacts and the estimated 
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value of g'^{A^). We also see a correlation between the fitting range and the resulting g'^{A^). But 
all values of g'^{A^) fall into the same ballpark and none of these fits can be done without such 
a positive contribution. To estimate the systematic uncertainty we have considered a large panel 
of fitting methods. All at more than four sigmas from zero except at /3 = 4.2 with the sliding 
window where they are only 2.5 sigmas above zero. Comparing the fitted {A^) using the tree level 
Wilson coefficient and that using the three loops one, we find that the latter is about 20 % below 
the former. 

The perturbative contribution to Zq, Z^'^^^ has a linear dependence in the bare lattice coupling: 



see Fig. [6jand Eq. (5.10), as expected from perturbation theory, but with a larger coefficient, even 
when the boosted coupling constant is used in perturbation theory. 

We also merge all three /3's after having subtracted the a^p^ term and rescaled the /3 = 4.05, 4.2 
to 3.9 using the ratios of Z^'^'^^{n). The overlap of the three data sets is rather good. The need 
of a non-perturbative contribution is also visible there. Both perturbative and non-perturbative 
contributions decrease with the momentum and are convex. This makes the separation difficult. 
Grossly speaking they share the decrease between 4 and 40 GeV^ in equal parts. 

We have converted our results for Z^^ and its perturbative part into the MS scheme. 

Combining all the results we find, using the three loop Wilson coefficient: 

g\^i'){A'),2cM = 2.01(11) (1°:?^) GeV^ = 10 GeV , 

ZP-*((2GeV)^^Le) = 0.773(3)(7) - 0.323(6) KB (^7Le- 1-5), 

jf{{2GeVr,glJ = 0.805(14) - 0.323(6) (^^2^) (^Le - 1-5) , 

P-*((2GeV)^^7Lc) = 0.750(3)(7) - 0.313(20) (^Le- 1-5) , 

^MSnon-perturbativc^^2GeV)^^/Lc) = 0. 781 (6) (21) - 0.313(20) ((?Lc " 1 -5) ; (6.1) 

where the systematic error is estimated from the scattering of the results in tables |3j [4], |5] and [7} 
We use the lattice spacings listed in table [T] 

Futhermore, table |8] also shows a nice agreement between the condensates for Nj = 2, although 
some systematics appears not to be under control for Nf = 0. This supports the interpretation of 
the oc 1/p^ contribution as being due to a condensate of the only dimension two operator in Landau 
gauge: A"^. Another confirmation comes from the validity of Martinelli-Sachrajda's criterium |10] . 
The accuracy on g'^{A^) is however limited due to several correlations in the fits. Further and 
more accurate checks of the consistency of g'^{A^) from other renormalisation constants will be 
very welcome. 
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A Appendix: The Wilson coefficients at 0{a^) 

The purpose of this appendix is to describe briefly the OPE analysis of the quark propagator 
renormahzation constant defined in Eq. (1.2 ) that leads us to Eq. (5.1 ), where the four-loops results 
in ref. [9] are exploited to derive the Wilson coefficients with the appropriate renormahzation 
prescription. This OPE analysis is analogous to the one performed in refs. [2||3||5||6] . The starting 
point is the OPE of the inverse of the quark propagator: 

C2Z /42\ 



zr\ii^) ' p' 4(iv2-i) 



ab 



(A.l) 



where only the leading term in ^ is kept, the quark mass being assumed to be negligible or to 
vanish. The cut-off regularization dependence is omitted for the bare quantities but that on 
the renormahzation momentum, /x, is explicitly written for the renormalized ones. In the RI'- 
MOM scheme we define Zf^' such that (^barc)"H/) = (p^) and such that, 

('^bare)~^(P^) ~ "^l^^afe Z^°^^ (p'^) . Then, the renormahzation momentum, /i^, taken to lie on the 
perturbative regime, one can apply Eq. (1.2) and obtain, 

+ C2Z, ' 



- Co.,(^,a(,)) +c..,(^,a(,)) ^^f^. (A.2) 

which implies a definition of cqz, and where C2Zg(p^,/i^) is the Wilson coefficient of g'^{A^). Al- 
though not yet specifying the renormahsation scheme of Z^'^'^^{fi'^), we know that 

cozJl,a(/^)) = l + 0(a2) , (A.3) 



while C2Zq is known up to 0{a'^) in the MS scheme j9j and, in particular, (1, a{^)) is given in 
eq (18) of that paper using = fi"^. Let us however keep in mind that 

C2zjl,a(/x)) = ^«(^) (l + 0(a(^))) = ^4^ (l + 0(a(^))) . (A.4) 
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Now, with the help of the appropriate renormahzation constants, one can also write Eq. (A. 2) in 
terms of bare quantities: 



(A.5) 



where Aj^ = . Then, as the /i-dependence of both l.h.s. and r.h.s. of Eq. (A.5) should 

match each other for any p, one can take the logarithmic derivative with respect to yU and infinite 
cut-off limit, term by term, on r.h.s. and obtains: 



7g(a(/i)) + 
-7A2(a(/i)) +7g(a(/i)) + 



d 



d log jj? 
d 



aiog;.^+^^"^^))^} 1-^2^, (^'"(/^) 








(A.6) 



where the /3- function, chosen to be in MS, is defined as 



d 



d log /x' 



1=0 ^ 



A-n 



i+2 



(A.7) 



and where 7g(a;(/i)) and 7^2 (Q;(/i)) are the anomalous dimensions for the fermion propagator and 
local operator A^, respectively, which are formally defined as 



7x(a(/i)) 



d 



d log fi' 



\ogZ 



X 



E 

i=0 



X 



Air 



(A.8) 



where X stands for q ot A^. The scheme for the anomalous dimension of is imposed through 
the renormahzation of the local operator A^, as was done in ref. [2] to obtain its leading logarithm 
contribution, and it is only known in the MS at the order 0{a'^) 
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Then, that is the only 

possible choice of scheme for 7^2 in eqs. (A.6). Concerning Z^'^'^^, its scheme is determined by the 



renormahzation prescription for the non-perturbative propagator in the left hand-side of Eq. ( A.l ) 



Both MS and RI'-MOM are possible. As we aim to obtain a non-perturbative formula to be 
confronted to the lattice estimate of the RI'-MOM quark renormahzation constant, it is convenient 



also to prescribe the RI'-MOM scheme for Z^'^^^. Thus, eqs. (A.6) must be re-written as 



7fl 



d 



d 



d log yU^ da 
a log oa 



In q 



Ri' r 



In c 



,WMS 

2Za 



0, (A.9) 
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where 



.Ri' P 



^pertRI'^p2) 



(A.IO) 



and c^^^ is in the WMS scheme^ exphcitly defined by the second equation of (A. 9), after the 
RI'-MOM prescription for Z^'^^^, that of MS for and by the choice of a boundary condition, 

„WMS, ' ' 



l,a{q)). Then, from Eq. (A. 2) one obtains 



.Ri' r 

0Z„ I ^2 



.RI' 

-OZn 



(A.ll) 



and, in practice, both eqs. (A. 9) can be combined to give the following differential equation, 

d 



d log /i^ 



d_ 
da 



.WMS / P 



2 ) 



.RI' 







(A.12) 



that can be solved to provide us with the ratio of Wilson coefficients, C2z,/coZg5 required to 
implement Eq. (|AAj_j). For the purpose of the best comparison with the results from the analysis 
performed in ref. we applied 



c™(l,a(p)) ^ c^|jl,a(p)) 



MS 

1 



(A.13) 



where (1, is taken from eq. (18) of ref.. 9 using if = /i^, as a boundary condition which 

is equivalent to the one applied in the analysis ofref. [g]. 

On the other hand, if we take Z^'^''* to be renormalized in MS, the equations in (A. 6) reads 



7r(«(/^)) + 



-7^^'(«(/i))+7r(«(/^))+ 



9 n/ / ^^ 

^ + a /i — 

a log jji-^ da 



In q 



MS / 1 



0Z„ 



, Q;(/i) 



d 



d log /i^ 



d_ 
da 



In c 



MS / ^ 



2Za 



2' 



Q;(/i) 



0, (A.14) 



where is the Wilson coefficient computed in ref. j9 , provided that the boundary condition, 
(1, a;(/i)), is taken again from eq. (18) of the same paper using = yU^. Then, we can 

^We define this scheme by imposing that the local operator of the Wilson expansion is renormalized in MS, 
while the expanded operator (the quark propagator, in our case) is in a MOM scheme. We called this "Wilson 
MS". 
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combine again Eqs. (A. 14) to obtain for c^^/c^^ the same equation Eq. (A. 12) that, with the 
same boundary condition, leads to: 



.WMS I P 



r,a(/i) 



"OZa 



r,a(/i) 



.MS 



r,a(/i) 



(A.15) 



On the other hand, we can also combine the second equation of (A. 9) with the second one of 



Eq. (A. 14) and obtain 



7f («(/^))-7r(«(/^))+ 



d 



d 



(91og/i^ ' ^ ^'~''da 



c 



WMS / P 



2Za 



2 ) 



MS 

^2Za 



(A.16) 



that, according to Eq. (A. 13), can be solved with the boundary condition ( 1 , o; (p) ) / ( 1 , a (p) ) 

1 and leaves us with a relation of WMS and MS Wilson coefficients wich allows Eq. (A. 11) to be 
re-written as 



c: 



X 



1 + 



MS I P 

2Za 



.WMS 

■^2Za 



MS,m2 



^RI' 



,2 5 



c 



MS / p 



(A.17) 



where, futhermore, is to be taken from ref. |91 and Cni' from ref. 



4(A^2 _ l)p2 

. Thus, we can use either 



13 



Eq. (A. 11) with the solution of Eq. (A. 9) or Eq. (A.17) with that of Eq. (A.16) to confront to the 
lattice estimates of Z^^ . Both expressions are equivalent. In the first case, one can proceed as 
was done in ref. [T] to solve Eq. (A. 9). To illustrate this first method, let us remind that Eq. (A. 9) 
can be solved at the leading logarithm by applying the following ansatz. 



.WMS / P 



327r 



.Ri' r 



a{p) 



(A.18) 



where we apply Eq. (A. 4) and the exponent a, required to satisfy Eq. (A. 9), should be 

105 - 8A^/ 



76 
/5o 



132 - 8A^ 



(A.19) 



/ 
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In the second case, to solve Eq. (A. 16), a similar ansatz extended to three-loops order can be 
applied. 



^WMS 













a{p) 



I 



1 + 



47r 



47r 



(A.20) 



where we use Eq. (A. 13) for the boundary condition. Then, by requiring that the ansatz Eq. (A.20) 
verifies Eq. (A. 16), the coefficients b and r^'s will be obtained in terms of those for the fermion 
propagator MS and RI'-MOM anomalous dimensions and for the MS /3-function. However, in this 
case. 



qMS qRI' 

7o -7o 

f3o 







(A.21) 



because the first-loop coefficient for the anomalous dimension is scheme independent (in the partic- 
ular Landau gauge, this scheme-independent first- loop coefficient is also zero for any scheme pL3]). 
Furthermore, as can be seen in the appendix C of ref. 
Landau gauge. Then, 



13], one is also left with 7^^^ = 7^ 



gRI' 



m 



ri 



Wo 







(A.22) 



and the Wilson coefficients for MS and RI'-MOM will thus differ only at the order O (a^), with 
the non-zero r^'s coefficients to be applied in Eq. (A.20) given by 



r2 



qMS qRV 

72 - 72 
_2/3o 

qUS qRV 

73 -73 



-25.4642 + 2.3333 A^^ , 



(A.23) 



/3: 



oMS qRV 

72 - 72 



-1489.9796 + 246.4424 Nf - 6.4609 



where the three and four-loop coefficient in MS and RI'-MOM for the fermion propagator anoma- 
lous dimension have been again obtained from ref. 13 . This leads, using Eqs. ( A.17[A.20 A.23) 
with A^^ _ 1 = 8 and Nj = 2, to our final formulae for the free-of-artefacts lattice determination 

of Zn 



^q ■ 



7Latt artefree/'„2 



/2 



(A.24) 



1 + 



.MS 



($'«(/^)) 1- 0.1317 0.5155 (^^2^ 



MS,At2 



^ c^l (^4, a(/i)) 1 - 0-1317 a'^ip) - 0.5155 a%p) 32 p'^ 
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In this last equation, we exploited the fact that the expression in parenthesis in Eqs. ( A.17[[Al24 ) 
does not vary with the renormalization momentum for the local operator A^, as can be inferred 
from Eq. ( A.12 ). Thus, once a gi ven m omentum, /x'^, is fixed for the renormalization of the fermion 
propagator in the l.h.s. of Eq. (A.l), the one appearing in Z^^"^^ in front of the expression in 
parenthesis, one is still left with the the freedom of choosing a renormalization momentum, /i^, 
which does not n eed to be the same, for the local operator inside the parenthesis. 

In Eq. (A. 24) the coefficients c^J^ and are known from perturbation theory, the former 
can be obtained from ref. 
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and the latter from ref. 



|9]. Two parameters are to be fitted: 
It is important to underline that 

the condensate is defined via the OPE, i.e. from Eqs. (AAWA^24-). Its precise definition depends 



non perturbative condensate (7 (/i) )_r. 



on the renormalisation scheme,the renormalisation scale, as well as the order in perturbation theory 
to which the coefficients Cqz, and C2z, are used. In Eq. (A. 17) and Eq. (A. 24) the renormalisation 
scheme for g^^ij) {A)r^^2 is MS and the scale is /i (10 GeV in our calculations). The coupling 
we use for the perturbative expansions of these coefficients, Cqz^ and C2Zq, is also chosen to be 
the MS one. These choices are kept all along the present paper. If we now wish to compare 
g^{lj) {A^)ji^^2 from the present calculation to that from another calculation, for example from 
the strong coupling constant [T], we must as far as possible use the same precise definition in both 
cases. However, its dependence on the scheme and on the order in perturbation theory is not so 
important: as seen in section |5] other systematic uncertainties are larger. 
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